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We study in detail the spectrum of heavy quarkonia with different orbital angular momentum along with their 
radial and gluonic excitations. Using an anisotropic formulation of Lattice QCD we achieved an unprecedented 
control over statistical errors and were able to study systematic errors such as lattice spacing artefacts, finite 
volume effects and relativistic corrections. First results on the spin structure in heavy hybrids are also presented. 



1. INTRODUCTION 

Heavy quarkonia are an interesting testing 
ground for our theoretical understanding of QCD. 
Many models have been developed to describe 
the vast amount of experimental data, but more 
accurate predictions from first principles will be 
needed for future experiments with heavy quark 
systems. However, the conventional lattice ap- 
proach becomes overly expensive if one tries to 
accomodate all the non-relativistic energy scales 
on a single grid: Mv 2 <C Mv <C M. Here v is the 
small velocity of the heavy quark with mass M. 
The non-relativistic nature of such problems has 
frequently been employed to formulate effective 
theories in which spatial and temporal directions 
are treated differently at the level of the quark 
action 00]. In particular the NRQCD approach 
has lead to very precise calculations of the low 
lying spectrum in heavy quarkonia || 

Additional problems arise if high energetic ex- 
citations are to be resolved on isotropic lattices. 
Very often the temporal discretisation is too 
coarse and the correlator of such heavy states can- 
not be measured accurately for long times. Glue- 
balls and hybrid states are among the prominent 
examples of such high-lying excitations and they 
receive much theoretical and experimental atten- 
tion as they are non-perturbative revelations of 
the gluon degrees of freedom in QCD. Early lat- 
tice studies have predicted such states starting at 



*A. AH Khan, S. Aoki, R. Burkhalter, S. Ejiri, 
M. Fukugita, S. Hashimoto, N. Ishizuka, Y. Iwasaki, 
K. Kanaya, T. Kaneko, Y. Kuramashi, T. M., K. Nagai, 
M. Okawa, H.P. Shanahan, A. Ukawa and T. Yoshie 



around 1.4 GeV, but those results were spoilt by 
large statistical uncertainties . This suggests 
that the inverse lattice spacing should at least be 
3 GeV or higher. More recently anisotropic lat- 
tices have been used to circumvent this problem 
in glueball calculations by giving the lattice a fine 
temporal resolution whilst maintaining a coarse 
discretisation in the spatial direction || . The suc- 
cess of this approach has triggered new efforts to 
measure hybrid potentials in an anisotropic gluon 
background and a comparative NRQCD analysis 
has demonstrated the validity of this adiabatic 
approximation for bbg hybrids. These attempts 
were reviewed in Q . 

In a previous study we reported on first quanti- 
tative results for charmonium and bottomonium 
hybrid states from anisotropic lattices p0| . Here 
we extend those methods to study also other exci- 
tations and the spin structure in heavy quarkonia 
more carefully. 

In Section 2 we present the details of our calcu- 
lation and results for the spin-averaged spectrum. 
In Section 3 we investigate the spin structure in 
heavy quarkonia and report on novel results for 
the splittings in heavy hybrids and D-states. 



2. SPIN-AVERAGED SPECTRUM 

In order to study excited states with small sta- 
tistical errors it is mandatory to have a fine res- 
olution in the temporal lattice direction, along 
which we measure the multi-exponential decay 
of meson correlators. To this end we employ an 
anisotropic and spatially coarse gluon action: 
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Here (j3, £) are two parameters, which determine 
the gauge coupling and the anisotropy of the lat- 
tice. Action (Q) is Symanzik-improved and in- 
volves plaquette terms, P M „, as well as rectan- 
gles, R^v It is designed to be accurate up to 
O(a*,o£), classically. To reduce the radiative 
corrections we invoked mean-field improvement 
and divided all spatial and temporal links by 
a 'tadpole' coefficient uq s and uot, respectively. 
We determined those coefficents self-consistently 
by measuring spatial and temporal plaquettes: 
uos = (Tr Py) 1 / 4 and u ot = (Tr P it ) 1 /*. With 
this prescription we expect only small deviations 
of £ from its tree-level value a s /a t . 

To describe the forward propagation of heavy 
quarks in the gluon background (An) we used the 
NRQCD approach introduced in H: 
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Here the NRQCD Hamiltonian, H, is designed 
to account for relativistic corrections and in- 
cludes spin-dependent operators up to 0(mv 6 ). 
The implementation details of Equation (Q) on 
anisotropic lattices are given in JllJ] . Since we are 
working with spatially coarse lattices it is cru- 
cial to improve all lattices derivatives and colour- 
electromagnetic fields in Equation (^|). Following 
the prescription of |jj we also achieved an accu- 
racy of O(o,g, a%) in the quark sector. 

At each value of the coupling, /3, we carefully 
tuned the heavy quark mass, M, so as to repro- 
duce the experimental ratios Mki n / (IP— IS) very 
accurately. Fortunately, the spin-independent 
quantities, such as IP — IS, are not very sen- 
sitive to the actual value of M, but the spin 
structure will show a strong dependence. For 
our calculation we chose several different values 
of the coupling, 0, which correspond to spatial 
lattice spacings between 0.15 fm and 0.47 fm. On 



even coarser lattices we cannot expect to control 
discretisation errors with our simple-minded ap- 
proach, while much finer lattices will violate the 
validity of NRQCD which requires a s M > 1. 

From the quark propagator, Gt, we construct 
meson correlators for bound states with spin 
S = (0, 1) and orbital angular momentum L = 
(0,1,2). Magnetic hybrid states are constructed 
from the colour-magnetic field coupled to a QQ- 
pair (Bi = [Aj, A/j]). For example, the spin- 
singlet operators read 

Q f Q , Q^iQ , Q^AjA k Q and B t Q . (3) 

Those simple operators can be further improved 
upon in order to optimise their overlap with the 
states of interest. Here we extract the excitation 
energies from multi-exponential fits to several dif- 
ferent correlators. 

Within the NRQCD approach one cannot 
extrapolate to the continuum limit and it is 
paramount to establish a scaling region for physi- 
cal quantities already at finite lattice spacing. In 
a previous study we found such scaling windows 
for the spin-averaged gluon excitations in both 
charmonium and bottomonium [ |Uj . These re- 
sults are particularly encouraging as they are in 
excellent agreement with calculations on isotropic 
lattices HQ], but with much smaller errors. In 
addition we were able to check our predictions 
against possible systematic errors such as finite 
volume effects. This is a natural concern since hy- 
brids are expected to be rather large owing to the 
flat potentials they are living in. In Figure |l| we 
show our results from lattices all larger than 1.2 
fm in extent, beyond which we could not resolve 
any volume dependence. For the Bottomonium 
hybrid we also have consistent results from two 
different anisotropics (£ = 3, 5) which confirms 
our initial assumption of small temporal lattice 
spacing artefacts. In fact, on our lattices the tad- 
pole coefficient uot deviates from its continuum 
value by only about 5% or less. It is also in- 
teresting to notice the presence of scaling viola- 
tions which are clearly visible if the lattices are 
too coarse for the physical system. 

We have now extended our analysis to study 
also higher radial excitations and D-states with 
L = 2. The spin-independent results are sum- 
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Figure 1. Scaling analysis for spin-averaged hy- 
brids. As we only measure excitation energies rel- 
ative to the ground state it is natural to present 
our results as the ratio Rb = (1£? — IS) /(IP — 
IS), which gives the normalized splitting of the 
magnetic hybrid excitation above the IS. 

marised in Figure |2[ The possibility to resolve 
all these excitations reliably should be considered 
the main success of anisotropic lattices. We are 
presently performing a finite volume analysis of 
the higher radial excitations. However, with the 
newly achieved accuracy we can also study spin- 
splittings in more detail. 

3. SPIN STRUCTURE 

Our inclusion of relativistic corrections to 
Equation (Q) is a significant improvement over 
previous NRQCD calculations of hybrid states, 
which were restricted to only leading order in 
the velocity expansion: (D(mv 2 ). At this level 
there are no spin-dependent operators and we 
have a strict degeneracy of all singlet and triplet 
states. Within the NRQCD framework, the spin- 
dependent operators appear first as higher order 
corrections to the Hamiltonan: 0(mv 4 ). This is 
in accordance with the experimental observation 
that spin-splittings in quarkonia are suppressed 
by v 2 compared to the spin-independent structure 
discussed in the previous section. Here we also 
include spin correction terms up to 0(mv 6 ), and 



Figure 2. Scaling analysis for excited Quarkonia. 
We plot the ratio Rx against the spatial lattice 
spacing for different states X = S,P,D and 3 B± + 
(exotic hybrid). 

study the breaking of the degeneracy. In partic- 
ular, we could directly observe the exotic hybrid, 
1 '", which is the state of greatest phenomeno- 
logical interest. 

Our results for fine structure and hyperfine 
splittings are shown in Figure || and ||, respec- 
tively. There are several interesting observations 
one can make. First of all, we find a notice- 
able reduction of the fine structure in D-states 
when compared to that of P-states. Similarly, 
the hyperfine splitting between spin-singlet and 
spin-triplet states is equally suppressed as the or- 
bital angular momentum is increased. This is 
in accordance with potential models which pre- 
dict a hyperfine splitting for only (L=0)-states 
since all other wavefunctions vanish at the ori- 
gin. Our data also indicates that the fine struc- 
ture in hybrid states is enlarged compared to the 
splittings in P-states, but we could not yet re- 
solve any splitting between the spin-triplet state 
( 3 B = 5 3 B 2 + 3 3 B 1 + 3 B Q ) and the singlet 1 B 1 . 

Finally one should also notice that the scaling 
behaviour of the spin structure is more involved 
than that of the spin-independent spectrum. This 
is not surprising since we have adopted a very 
simplicistic approach to determine all the coeffi- 
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Figure 3. Fine Structure in Charmonium and 
Bottomonium on three different lattices. All re- 
sults come from NRQCD calculations with accu- 
racy 0(mv e , a,g, a\ ). The conversion to dimen- 
sionful units is done by setting the IP — IS split- 
ting to its experimental value. 

cients in NRQCD with a single prescription (tree- 
level tadpole improvment). Namely the hyper- 
fine splitting 3 Si — 1 S'o does not scale on the lat- 
tices considered here. It is apparent that one 
needs a better improvement description to ac- 
count for lattice spacing artefacts in such UV- 
sensitive quantities. 

In conclusion, we find that coarse and 
anisotropic lattices are extremely useful for preci- 
sion measurements of higher excited states. This 
is due to an improved resolution in the tempo- 
ral direction and the possibility to generate large 
ensembles of gauge field configurations at small 
computational cost. It has lead to an unprece- 
dented control over statistical and systematic er- 
rors in lattice studies of heavy quarkonia. Af- 
ter the inclusion of relativistic corrections we are 
also sensitive to spin-spin and spin-orbit interac- 
tions and could observe a clear hierachy in the 
spin structure, depending on the orbital angular 
momentum. Spin splittings in hybrid states were 
found to be larger as a result of the gluon angular 
momentum to which the spin can couple. The re- 
maining and dominating systematic error for all 
our predicitions is an uncertainty in the scale as 



Figure 4. Hyperfine splitting between the aver- 
aged spin-triplet and spin-singlet states of differ- 
ent orbital angular momentum. 

the result of the quenched approximation. This is 
not yet controlled and we find a variation of 10- 
20%, depending on which experimental quantity 
is used to set the scale. 
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